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ABSTRACT 

Lyman a emission from high-redshift galaxies may be a powerful probe of the ioniza- 
tion history of the IGM at z > 6: the observed Lyman a emission line is sensitive to 
the neutral fraction of IGM hydrogen in the range 0.1-1. We present calculations of 
observed Lyman a emission lines from z > 6 galaxies, illustrating the effect of varying 
the many free parameters associated with the emitting galaxy, its halo, and the IGM 
around the galaxy. In particular, we use a dynamic model of the IGM that includes 
the effect of IGM infall toward the emitting galaxy. Galactic winds may play a crucial 
role in determining observed Lyman a line fluxes. We compare our model predictions 
with observations of two z = 6.5 galaxies and conclude that, if galactic winds are al- 
lowed for, existing observations place no constraint on the neutral fraction of the IGM 
at z = 6.5. Future space-based observations will constrain the importance of galactic 
winds; if winds are unimportant for the observed z = 6.5 galaxies, our models suggest 
that the IGM neutral fraction a± z — 6.5 is < 0.1. 

Key words: line: profiles - intergalactic medium - galaxies: high-redshift - diffuse 
radiation - cosmology: theory. 



1 INTRODUCTION 

Strong Lyman a emission is present in many distant galax- 
ies, and is an important signpost for discovering high- 
redshift galaxies and, especially, measuring redshifts for 
them. The presence of strong Lyman a emission has been 
crucial in confirming the highest redshift galaxie s known 
at g > 5 (e.g., FHu. McMahon. & Cowic 1999; Hu et al 



2002a b; Aiiki ct al. 2003; Kodai ra et al.l2003ll:lRhoads et al 
i200a') . See .Taniguchi et al. (.2003) for a review. It also plays 
an important role in identifying the redshift of some galaxies 
at 2 > 2 because its strength, even in some apparently ex- 
tremely dusty galaxies, allows for redshift d eterminations 
when little continuum ligh t is visible (e.g., ISteidel et all 
l2003l:IChapman et ahlliooa) . 

In this paper we concentrate on Lyman a emission pow- 
ered by star formation. Hot stars emit photons capable of 
ionizing hydrogen. If the ionized gas density is sufficient, 
hydrogen will recombine on relevant timescales. The recom- 
bination cascade produces a rich spect rum of lines, but usu- 
ally includes a Lyman a photon (see IOsterbrockl|l989l . for 
a review). Because most recombinations produce Lyman a 
photons, it is predicted to be the strongest line emitted from 
recombination emission. AGN activity can also provide ion- 
izing photons that ultimately generate a strong Lyman a 
line. Additional sources of L yman a photons are atomic- 
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1200(1 iFardal et alll200lD and fluorescence of gas clouds il- 
luminated by a stro ng, non-local source of ionization (e.g., 
iReuland et"alll2003l) . 

The predicted strength of the Lyman a line, and the 
assumed prevalence of star formation at early times (in- 
ferred from the abundance of old stars in the local uni- 
verse) led to an early prediction of an abundant population 
of high-re dshift galaxies with very s trong Lyman a emis- 
sion lines JPartridee fc Peebles! IT967I) . Though subsequent 
revisions to the theory of galaxy formation, primarily that 
it was a hierarchical process, explained the relative lack of 
Lyman a detections com pared to those early predictions 
iHaiman fc SpaanJl999^ , searching for galaxies based solely 
on strong Lyman a emission has been validated recently as 
an effective technique for identifying galaxies up to z = 6.5. 

Lyman a is a resonant transition with a large cross- 
section. As a consequence, small quantities of neutral hy- 
drogen scatter away Lyman a photons in direction and fre- 
quency. In particular, if the intergalactic medium (IGM) 
is not highly ionized, as is expected for at least some of 
the time between recombination (z ~ 1100) and z ~ 6 
jFan et alll200l iKoeut et alJliooa) . then the neutral hy- 
drogen in the IGM can easily scatter the Lyman a line from 
galaxies. The Lyman a photons are not destroyed, but the 
scatteri ng process diffuses the line in areal coverage and fre- 
quency jLoeb fc Rvbickll999ll . and the resulting low surface 
brightness emission is currently unobservable. 

A simple consideration of the IGM shows that the blue 
side of the emergent Lyman a line would be scattered away 
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by a neutral IGM, and only the red side of the line would 
be observed. This follows from considering that most of the 
scattering by neutral hydrogen comes at wavelengths near 
the Lyman a transition in the rest-frame of the gas. How- 
ever, for a completely neutral IGM at high redshifts, the Ly- 
man a cross-section far from resonance becomes important, 
since the natural line profile has 'damping wings' that fall 
off only like z/^^. Thus even if the IGM has no component at 
the velocity of corresponding to resonant scattering of a par- 
ticular Lyman a photon, that photon may still be scattered. 
This observation led to the conclusion that observations 
of Lyman a emission lines (or their absence) at high red- 
shifts may probe the ionization state of the IGM in a regime 
poorly test ed by the Gunn-Peters o n trough measurements, 
xm ~ 1 I'Miralda-Escud el ll99S|: [ M iralda-Escude fc R.eesI 
ll998l: |Haiman 2002; Barka na fc Loebi 2003b: Con 2003b). 

In particular, since the WMAP satellite discovered a 
high optical depth to the last scattering surface (iKoeut et alJ 
l2n03r . the ionization history of the universe at 2 > 6 
is of great interest. Though the WMAP results indicate 
that the universe was mostly ionized for a substantial his- 
tory of the universe at z > 6, it does not constrain 
the exact ionization level or the history in detail. There 
are now many predictions for the ionization history at 
z > Q based on semi-an alytic and numerical modelling of 
early star formation (e .g., Ccn 2003a; Wvithe & Locb 2003; 
ICiardi. Ferrara. fc White ,2003 : Haiman fc Holder 2003) . If 
there is a population of strong Lyman a emitters at the red- 
shifts in question, 6 < z < 20, and if the effect of IGM scat- 
tering on their observed properties can be calculated, then 
future surveys for Lyman a emission may provide valuable 
information on the reionization history at z > 6. 

Here we use a more sophisticated treatment of the dy- 
namics of the IGM around galaxies with strong intrinsic Ly- 
man Q emission to investigate the effect of the IGM on the 
observed line. The IGM model we use will be presented in 
more detail separately, but i s summar i zed h ere and nearly 
identical to that described in iBarkanal ||20o3)- The key fea- 
ture is that infall of the IG M toward the galaxies is included. 
iBarkana fc Loebl ll2003al) have showed that such an IGM 
model makes a prediction for an absorption feature in the 
spectra of high-redshift QSO Lyman a lines. Here we con- 
sider normal galaxies, where the effect of an infalling IGM 
is even stronger. 

Throughout this paper we assume a cosmology based 
on the recent res ults from WMAP and other work 
iSpergel et alJ lioOSl and references therein): 57m = 0.3, 

= 0.7, Ho = 70 kms"\ and f2b = 0.0469. In general 
we quote distances in physical Mpc, denoted by 'pMpc' 

This paper is organized as follows. In Section |5| we 
present the Lyman a scattering cross-section. Section |3 de- 
scribes the model for the intrinsic Lyman a line profile and 
the observed profile modified by IGM scattering. Section 2] 
presents Lyman a line profiles as a function of the (many) in- 
put parameters that influence the observed line profile. We 
introduce the possible role of galactic winds in Section |S| 
Section 15] discusses our Lyman a flux predictions with com- 
parison to current observations, and also reviews model as- 
sumptions. Section summarizes. 



2 LYMAN Q SCATTERING 

A Lyman a photon is emitted when a hydrogen atom makes 
a transition from the n = 2 level to the n — 1 ground-state. 
The emitted photon has an energy of 10.199 eV and a wave- 
length of Aq = 1215.67 A. The transition is the strongest 
spontaneous transition of a hydrogen atom, with an Ein- 
stein A value of A21 = 6.265 x 10* s~^. Correspondingly, 
the absorption cross-section of the transition is relatively 
large as well. 

A Lyman a photon travelling through a region of neu- 
tral hydrogen will be repeatedly absorbed and re-emitted. 
This process by itself does not destroy Lyman a photons, 
but if photons are scattered over a sufficiently large region 
such that an observation spatially resolves the emission, then 
the Lyman a line wil l be difficult to observe in practice 
jLoeb fc Rvbickilll999ll . In this work we assume that Ly- 
man a photons scattered in the IGM, that is, outside of the 
virial radius of the emitting galaxy, are missed by observa- 
tion. For a 2 > 6 galaxy, the virial radius is larger than 1" 
(the size of a typical spectroscopic slit) for any galaxy more 
massive than ~ 3 x 10® M0. 

So under the assumption that photons emitted in the 
Lyman a emission line of a galaxy are lost from subsequent 
observation if they are scattered by the IGM, we constructed 
a model of the IGM (see Section I3.2.H and computed the 
absorption due to neutral hydrogen. The natural absorption 
cross-section ctn is lPeebles..l993.'l 



(Tn(«^) = 



8n ^■K^{v-v^Y + {AlJ^){vlv^Y' 



(1) 



where v = c/A and Va = c/\a. 

The total absorption cross-section due to a parcel of 
IGM gas is determined by the kinetic properties of the gas, 
both bulk and thermal, as well as the natural cross-section. 
Bulk motions just introduce a frequency shift, but thermal 
motions broaden the profile, modelled as the convolution of 
the natural cross-section with a Maxwellian velocity distri- 
bution, 



M{v) 



\2nkT) 



1/2 



exp 



muv 
' 2kT 



(2) 



(3) 



where mn is the mass of the hydrogen atom and T is the 
temperature of the IGM gas (modelled in the next section) . 
For the usual Lorentzian approximation of (tn, the convolved 
cross-section would be a Voigt profile, hence we label our 
convolved cross-section ay. Figure^shows ay{X), assuming 
a temperature of T = 10** K for the scattering gas. The 
profile has a Doppler core with a half- width of ~ 10 kms~^ 
and power- law tails called 'damping wings.' If the IGM at 
the redshift of interest is or has previously been reionized, 
then T = 10* K is probably a good approximation. If the 
IGM is still thermally pristine, then the temperature may 
be as low as ~ 20 K, resulting in a narrower Doppler core 
of the scattering proflle. In practice this would have a small 
effect on any near-future observation. 
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Figure 1. Cross-section for Lyman a absorption, as a function 
of wavelength difference from resonance in the frame of the ab- 
sorbing atom. The top axis shows the velocity (in the frame of 
the absorbing atom) where the Lyman a resonance corresponds 
to the wavelength on the bottom axis. 



3 LYMAN a LINE PROFILES 

We will refer to two types of Lyman a lines from galax- 
ies. The intrinsic line is the line produced in the galaxy, 
specifically, the line that reaches the virial radius of the halo 
containing the emitting galaxy. The observed line is the 
line after propagation though the IGM to z = 0. 



3.1 Intrinsic Lyman a emission line 

A Lyman a emission line from a high-redshift galaxy is typ- 
ically a consequence of hydrogen recombinations in nebular 
gas. The ionization of the H ll region may be maintained 
either by hot stars or by gas accretion onto a black hole. 
Though luminou s QSOs have been discovered at z > 6 
jFan et ahlliboj . and references therein), deep x-ray data 
indicates that most galaxies at those re dshifts do not con- 
tain powerful AGNs llBarger et aljl20o3) . Thus in this work 
we will assume recent (< 10 Myr ago) star formation powers 
Lyman a emission lines. 

Type O stars emit hydrogen-ionizing photons.^ Some 
fraction, /osc, of these photons may escape the galaxy with- 
out absorption. The rest of the photons are absorbed by ei- 
ther neutral hydrogen or dust within the galaxy. An ionizing 
photon absorbed by a hydrogen atom initiates a cascade of 
recombination emission, ultimately resulting in a Lyman a 
photon about 2/3 of t he time, and two -photon emission the 
other 1 /3 of the time iOsterbrockll989l) . A Lyman a photon 



^ Approximately 1 Gyr after the onset of star formation hot white 
dwarfs form, and the se may emit substantial ionizing radiation 
JCharlot fc Fal l''2000^: however, for z > 5.5 the universe is not 
old enough to have formed these stars. 



subsequently scatters through absorption and reemission un- 
til the optical depth of the galaxy at its wavelength is < 1; 
the optical depth is determined both by the position and 
direction of the photon, and also its wavelength. If dust is 
present, it may absorb photons at any time in this process. 

The radiative transfer problem for Lyman a depends 
sensitively on the geometry of neutral hydrogen and dust 

and also the dynamics of 
For simplicity we will as- 
sume that the intrinsic Lyman a emission line is centered at 
the systemic redshift of the emitting galaxy, with a Doppler 
profile shape described by a characteristic velocity, vu, 



within the gal axv (Neufeld 1991.) . 
the hydrogen jKunth et aljil99ai . 



exp 



(Ai^d) 



where 0"^ is the line profile per unit frequency and 



C 



(4) 



(5) 



In Sectionl^we address other possible intrinsic line profiles. 

In this simplified model, two parameters characterize an 
intrinsic Lyman a emission line: its luminosity, La, and vu- 
La depends on the rate of production of ionizing photons, Q, 
the ionizing photon escape fraction, /osc, and dust absorp- 
tion. The production of ionizing photons in turn depends on 
the star-formation rate and the metallicity and initial mass 
function (IMF) of the stars formed. None of these quantities 
may be reliably predicted from first principles for a galaxy 
at high redshift. Therefore, any use of the observed Lyman a 
line of a high-redshift galaxy to deduce the radiative trans- 
fer properties of the IGM will be complicated, unless these 
quantities can be measured independently of the Lyman a 
line. 

The estimation of t^D is simpler. The likely minimum 
velocity scale for motion in a galaxy is the circular velocity of 
the parent halo, v^; this velocity may describe the gas motion 
is the star-formation is mediated by a violent process, such 
as the merging of two equal-mass galaxies. If star-formation 
is a more quiescent process in the galaxy, it may occur in a 
disk whose peak velocity is, roughly , between V c and 2vc , for 
realistic halo an d disk properties tMo. Mao, fc White.il99^ : 
ICole et aLlboodl . Thus we assume Vc vo '^^ 2vc- 

3.2 Observed Lyman a emission line 

The optical depth for photons to be absorbed in the IGM is 



r(^) = 



1 + 



xm{r)nii{r)aY{u')dr, 
v{r) 



(6) 



(7) 



where -Robs is the physical distance to the observer (though 
there is almost no contribution from r ^ iivir because the 
Hubble flow eventually redshifts the IGM Lyman a reso- 
nance far from the emitted Lyman a line [in the frame of 
the galaxy]). We describe our model for nnir), v{r) and 
xm{r) in Section f3.2.1l 

The observed Lyman a emission line, in the rest-frame 
of the emitting galaxy, is then 

't'l^^'i'^) = 0'"(J^) exp [-r(i/)] . (8) 
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3.2.1 IGM model 

3.2.1.1 Density and velocity We will present a full de- 
scription of our model for the density and velocity of IGM 
material in another paper (Santos & Adelberger, in prepara- 
tion). The model is very similar to one described in lBarkanal 
l|2002l) . Here we summarize our model as it applies to this 
paper. 

The evolution of the density and velocity of the IGM fol- 
low from the small, linear, initial perturbations imprinted on 
the universe at very high redshift. We model the IGM start- 
ing from the statistical description of linear perturbations 
by applying a simple description of the non-linear evolution 
of the perturbations. 

Specifically, we start by considering a Lyman a emit- 
ting galaxy in a halo of total mass M at redshift z. We then 
constrain the initial linear overdensity on the scale of M 
to produce the halo at redshift z. Starting from that initial 
linear overdensity constrai nt, we use the excu rsion set for- 
malism t.Bond et al. 1991: .Lacev fc Cole|[^993^ to compute 
the mean'^ initial linear overdensities, averaged in spheres 
centered on the location of the halo, on scales larger than 
the scale of M. The mean initial linear overdensity decreases 
monotonically with increasing scale; this is not externally 
imposed. Our approach also implies that the halo is not sub - 
structure of a larger collapsed halo (see also lBarkanall20o3) . 

Once we have the initial linear overdensities of spheres 
enclosing the halo, we evo lve those overdensities usin g the 
spherical top-hat model jPartridee fc Peebles! ^^). Be- 
cause the mean initial linear overdensity decreases mono- 
tonically with increasing scale, there is no shell crossing, 
and the evolution of each IGM matter shell is independent 
of the evolution of the other shells. 

Our dynamical evolution treats all matter equally, and 
we assume the baryons trace the dark matter in these shells. 
This description would be inappropriate if, e.g., the Ly- 
man a emitting ga laxy blows a stro ng wind into its sur- 
rounding IGM (see lAdelbereer et ali .200.3. for evidence of 
this at 2 = 3), and we will return to this in Section |3 

A realistic distribution of initial linear overdensities 
around a halo, in contrast to the mean initial linear over- 
densities averaged over many halos, would not in general be 
spherically symmetric or even monotonic with radius within 
spherically-averaged shells. In addition to collapsing toward 
the galaxy, there would be structures collapsed and collaps- 
ing within the 'IGM.' The most straightforward way to as- 
sess the effect of these complicated dynamics would be with 
numerical simulation. Here we appeal to analytic results to 
demonstrate the effect of structure with the collapsing IGM 
will not severely affect our conclusions. 

Scannaoicco & Bar kana (20o3) computed the bivariate 
probability distributions of two halos forming at different 
positions. Their results suggest that outside the virial ra- 
dius of a halo, the typical enhancement in the number den- 
sity of halos is, at maximum, only a factor of a few, for any 
mass. Thus we conclude that the fraction of matter collapsed 
into halos (above a given mass) in the IGM surrounding our 
galaxy of interest is at most a few times the universal col- 
lapse fraction. At z = 6 up to half the mass of the IGM may 
be in the form of collapsed halos with virial temperatures 



above 10* K, that is, halos massive enough to collisionally 
ionize their hydrogen. That fraction decreases with increas- 
ing radius. The enhancement is also a decreasing function 
of decreasing mass. Thus we conservatively expect that our 
simple picture gives at most a factor of two overestimate 
of the optical depth, averaged over a reasonable sample of 
halos. At large optical depths, this error doesn't affect the 
observed Lyman a line properties, and at low optical depths 
we may underestimate the observed line by up to that same 
factor of two. 

The practical output of our model is, given a halo mass 
AI at redshift z, the radius, velocity, and mass of shells of 
IGM matter surrounding the galaxy from its virial radius to 
a large radius. 

3.2.1.2 Ionization state To calculate the scattering of 
a Lyman a emission line by the IGM, we need to know the 
ionization state of the hydrogen, in addition to its density 
and velocity profile. We solve the ionization balance of the 
IGM gas assuming two possible contributions, direct ioniza- 
tion by the galaxy emitting the Lyman a line and, if the 
universe is ionized, a mean ionizing background. 

After calculating the density distribution around a halo 
using the prescription described in the previous section, we 
solve the ionization balance as a function of distance from 
the galaxy. In general, the recombination rate of ionized hy- 
drogen in the IGM is relatively small. As a consequence, in 
the pre-reionized universe H ll regions around galaxies are 
not in equilibrium, but began growing when star-formation 
turned on and expand at the rate that ionizing photons can 
ionize more neutral IGM (rather than balancing recombi- 
nations within the ionized region). The Hll region is still 
expected to have a relatively sharp boundary, as the mean 
free path of ionizing photons in the neutral IGM is small. 
After reionization, all of the IGM is highly ionized by the 
mean ionizing background, but the direct ionizing flux of 
the galaxy provides an additional proximity effect on the 
ionization balance immediately around it. 

We divide the IGM around into one or two regions: ei- 
ther the universe is reionized, in which case all of the IGM is 
ionized; or the universe is not reionized, in which case there 
is an Hll region in the IGM immediately surrounding the 
galaxy, and the IGM is fully neutral outside of that region. 
In that case the radius Rs of the H ii region is calculated by 
solving 

'■''^ 4TTr^nn{r) 



n_ Q-A(r) 



-dr = tsf 



nH(r-) 



[Sir) + 1] pufi{l + zf 



— d/y, 
hv 



A(r) = / 47rr nH(f)aB(T)df, 



(9) 
(10) 
(11) 
(12) 



where _Rvir is the virial radius of the galaxy's halo, tsF is the 
age of the star-formation activity'', 5{r) is the overdensity 
with respect to the mean IGM density (calculated with the 



i.e., averaged over a large ensemble of identical halos 



We model star formation with a constant rate for an age tgp. 
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model in the previous section), pH,o is the comoving (i.e., 
z = 0) density of hydrogen, z is the redshift of the galaxy, 
hi^o is the ionization energy of hydrogen, is the specific 
luminosity of the galaxy, and qbCT) is the hydrogen Case 
B recombination coefficient. Since the Hll region is almost 
completely ionized, here we assumed the electron physical 
density, rie, is rie = wh, and ignored the contribution to ric 
from helium. We have not included a clumping factor mod- 
ification to the recombination rate; however, we do account 
for the mean overdensity of the IGM matter. Additional 
clumping of this material would decrease the size of the true 
Hll region compared to our calculation. 

The ionization rate for a hydrogen atom within an ion- 
ized region a distance r from the galaxy is 



r(r) 



hu 



(13) 



where Ji, is the mean specific intensity of the intergalactic 
background radiation. If the universe is not reionized, then 
we set Ji, = 0. The hydrogen neutral fraction at any radius, 
XHi{r), is solved for by setting r(r) equal to the recombina- 
tion rate at that radius, 



XHi(r) = l + |-i(x^ + 4x)^/' 



r(r) 



nH{r)aA{T) ' 



(14) 



(15) 



where aA{T) is the hydrogen Case A (total) recombination 
coefficient. Additional clumping of the IGM would decrease 
X linearly, increasing xm- 

We assume that when IGM gas falls across the virial 
radius it shocks and ionizes completely, and has no impact 
on the Lyman a line profile. 



4 LYMAN a LINE RESULTS 

As described in Section [3.11 the strength and shape of the 
intrinsic Lyman a emission line depends on many properties 
of the source galaxy. Additionally, the properties of the ob- 
served line are strongly modified by scattering in the IGM, 
which depends on the properties of the galaxy, its halo, and 
the IGM model. To illustrate these dependencies, we will 
adopt a fiducial model for a Lyman a emitting galaxy, then 
show the dependence of the intrinsic and observed Lyman a 
line on the properties of that galaxy and the IGM. 



4.1 Fiducial galaxy/halo model 

Our fiducial Lyman a emitter model is a galaxy at rest in 
the center of a 10 M© halo at z = 6.5. 'Fiducial galaxy' 
will refer to this total system, galaxy plus halo. We assign 
the fiducial galaxy a constant star-formation rate (SFR) of 
10 M0 yr~^ for an age of tsF = 10** yr. The virial radius of 
the halo is i?vir ~ 19.8 kpc, and the circular velocity at the 
virial radius is — 148 kms~^. The halo dynamical time is 
then 1.3 x 10* yr. Given our cosmology, and assuming the 
universal baryon to dark matter ratio holds in the halo, the 
halo contains 1.5 x 10^° Mq of baryons. Thus at 10 M© yr~^ 
the gas supply could last up to 1.5 x 10^ yr, about twice the 
age of the universe at that redshift. The chosen SFR also 



implies that 9 per cent of the baryons are converted into 
stars per halo dynamical time. 

We use the following prescription to convert the fiducial 
galaxy properties into the emitted Lyman a line. We convert 
the SFR to a hydrogen-ionizing photon luminosity using 



= 3.5 X lO'^^ s"^ 



SFR 

M© yr- 



(16) 



appropriate for a Salpeter IMF from 1 to 100 M© with 
l/20t h solar metallicity jLeitherer et al.l Il999l : ISchaereJ 
120031) . We assume that 2/3 of the ionizing phot ons that do 
not e scape are converted into Lyman a photons llOsterbrockl 
Il989h. so that 



(17) 



La = -{l- fcsc)QhVa- 

Note that this assumes there is no absorption by dust; we 
return to the issue of dust in Section r6.3l The other Lyman a 
line parameter that needs to be specified is the line width, 
characterized by «d = fv'^'c- For our fiducial model, we set 
U = 1.5. 

Additionally we need to model the effect of photoion- 
ization of the galaxy on the surrounding IGM. We assume 
the escaping photon luminosity of ionizing photons is simply 



CSC J CSC 



(18) 



Our baseline choice for the escape fraction is /esc = 0.1, 
giving 

SFR 



La = 3.4 X 10"*^ erg s"^ 



M© yr- 



(19) 



4.2 IGM model dependence 

Now we illustrate the IGM scattering of the Lyman a line 
emitted by our fiducial galaxy, highlighting the importance 
of using our dynamic IGM model. At 2 = 6.5, the ionization 
state of the IGM is unknown. Limits from QSO spectra at 
redshifts approaching 6.5 indicate that xhi > 10~^ averaged 
over the volume of th e IGM, and xbj > 10^^ averaged over 
the mass of the IGM llFan et alJl2002ll . Simulations of reion- 
ization suggest that the z ^ 6 transit ion from xm — 1 to 
Xm — was very sudden (|Gnedin"2000). Combined with the 
QSO observations, these suggest that xm could have almost 
any value ai z = 6.5. The WMAP results constrain the IGM 
ionization fraction to be ~ 1 over a substantial fraction of 
the evolution of the universe at 2; > 6, but do not provide 
detailed constraints on models of xm{z). 

We start by assuming our fiduci al galaxy (z = 6.5) 
is embedded in a fully neutral IGM. iHaimanl i2002h ar- 
gued that the detection of reasonably strong Lyman a emis- 
sion from 2 — 6.5 galaxies did not rule o ut a fully neutral 
IGM, as had been previously suggested iIHu et a lJ|2002a|: 
iRhoads et al.ll2'003^ . However, the calculations of lllaimanl 
(200'3) did not account for the infall profile of the IGM 
around a 2 = 6.5 galaxy. 

First we illustrate the difference between our model of 
the IGM, and a smooth, Hubble flow model. We call our 
model the dynamic IGM model because we allow the 
IGM to evolve away from the comoving solution. We call the 
mean-density, comoving model the simple IGM model. 
The simplest difference between the models is the density 
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Figure 2. IGM density as a function of radius from the fiducial 
galaxy. The left edges of the plots are at the virial radius, (a) 
The total IGM mass density and density contrast as a function 
of radius in physical Mpc. (b) The cumulative IGM hydrogen 
column density to radius r. In both plots the solid curve is for 
the dynamic IGM model and the dashed curve is for the simple 
model. 
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Figure 3. Velocity of IGM gas as a function of radius from the 
fiducial galaxy. The solid curve is the dynamic IGM model; the 
dashed curve is the Hubble flow relation of the simple model. 
The dotted curve is the ratio of the dynamic model to the simple 
model, in units marked on the right axis. 



of the IGM as a function of distance from the center of the 
halo. Figure shows the density and density contrast of 
the IGM outside the halo virial radius for our dynamic IGM 
model, and a simple IGM. Within about 10 virial radii, or 0.2 
physical Mpc (pMpc) , the density enhancement is between 2 
and 10 times the mean IGM density. Integrating the density 
out from the virial radius gives the curves shown in Fig.|2j3, 
after converting total density to hydrogen density. Though 
Fig. |5| shows that there is a substantial difference in the 
density profile at small radius, the increased density as a 
function of radius is not the most important effect of the 
dynamic IGM model. 

The most important difference between the IGM mod- 
els, with respect to scattering of Lyman a line photons, is 
the distribution of density as a function of velocity. Figure|3 
shows the relationship between radius and velocity, for the 
dynamic and simple IGM models. With this relation we re- 
plot the densities of Fig. |5| against velocity in Fig.|l| The lo- 
cal Hubble parameter is H{z — 6.5) — 790 kms~^ pMpc^^. 
For the IGM at positive velocities with respect to the galaxy, 
that is, the blue side of the line, there is little difference in 
density between the dynamic and simple models. However, 
while the simple model has no gas at negative velocities, cor- 
responding to the red side of the line, the dynamic model 
has reached a column density of hydrogen of lO^'' cm~^ by 
the time the velocity has reached 0. The absorption profile 
(Fig.0 has a cross-section of greater than 5 x 10~^* cm^ over 
the central 10 kms~^. Thus the average optical depth within 
about 100 kms~^ of the red side of the line is greater than 
5 X 10*xhi. Immediately we can see that for xm ~ 10~*, a 
typically expected value for the IGM, optical depths of order 
unity and greater may be expected. This optical depth on 
the red side of the observed line due to scattering by highly 
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Figure 4. IGM density as a function of velocity in the frame of 
the fiducial galaxy, (a) The total IGM mass density and density 
contrast as a function of velocity, (b) The cumulative IGM hy- 
drogen column density to velocity v. In both plots the solid curve 
is for the dynamic IGM model and the dashed curve is for the 
simple model. 

(but not completely) ionized gas gives rise to a line profile 
distinctly different from a simple IGM model. 

Figure |S| presents a schematic illustration of the differ- 
ence between the dynamic and simple IGM models. As de- 
scribed above, we see that the dynamic model has gas with 
negative velocity, in contrast to the strictly positive velocity 
IGM in the simple model. We also see that the Stromgren 
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Figure 5. Schematic diagram of location and velocity of the IGM 
around the fiducial galaxy, (a) Dynamic IGM model, (b) Simple 
model. The galaxy is marked by a 'G' at right. The curve im- 
mediately surrounding it shows the virial radius of the halo. The 
next curve, at ijg, shows the Stromgren sphere. The velocity axis 
marks the velocities (in kms~^) corresponding to these radii. In 
(a) -Rstall) the radius where IGM gas has zero velocity with re- 
spect to the galaxy, is also marked. The dashed and dotted lines 
key the physical location of the IGM material to the contribution 
it makes to the scattering optical depth, shown in Fig. |£] 



radius is both larger and extends to larger IGM velocities. In 
Fig.|^we have keyed diflterent physical regions with different 
line types; Fig. El shows the effect each of these regions has 
on the observed line profile. 

In Fig. El we show the optical depth as a function 
of observed wavelength, subdivided by the velocity of the 
IGM gas responsible for the optical depth, for the fiducial 
galaxy with the dynamic and simple IGM models. We la- 
bel the velocity of the IGM at radius r from the galaxy 
v{r); v{Rs) = 222 kms^^ for the dynamic model, and 
v{Rs) = 373 kms~^ for the simple model. When examining 
the observed line profile at a particular wavelength A, it is 
useful to define the velocity in the frame of the galaxy where 
the Lyman a resonance corresponds to that wavelength. 



1 



A 



il + z)K 



(20) 



The total optical depth has three prominent regions. 
At Va > v{Rs), the optical depth is dominated by Doppler 
core scattering from neutral IGM at ii ~ Hq,. At v{Rvii) < 
Va < v{Rs) the optical depth is still due to Doppler core 
scattering, but the IGM ai v ~ Va is mostly ionized. The 
third region is at Va < w(-Rvir), where the optical depth is 
due to IGM gas at v > v{Rs) (the same as the first region), 
but now the scattering results from the damping wing of the 
scattering cross-section. Figure |S1 illustrates that the central 
part of the observed line profile depends on the ionization 
state of the IGM gas in the Stromgren sphere, and the red 
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Figure 6. Optical depth due to IGM scattering as a function of 
observed wavelength, (a) Dynamic IGM model, (b) Simple model. 
The solid curves are the total optical depth due to all IGM H I. 
The short-dashed curves are the optical depth due to gas with 
negative velocity (in the frame of the galaxy). The dotted curves 
are the optical depth due to gas with positive velocity inside the 
Stromgren sphere. The long-dashed curve is the optical depth due 
to gas outside the Stromgren sphere. 



wing of the line profile depends on the size of the Stromgren 
sphere and the IGM outside of it. 

Figure|7|quantitatively compares the IGM density com- 
puted for both the dynamic and simple IGM models. In the 
dynamic model, Rs = 0.38 pMpc, versus Rs = 0.47 pMpc 
for the simple model. The larger Stromgren sphere in the 
simple model is a consequence of the lower IGM density. 

Figure |S1 shows the observed Lyman a line profiles per 
unit wavelength, (fix, for the dynamic (solid curve) and sim- 
ple (dashed curve) IGM models, compared to the intrinsic 
line profile (dotted curve) of the fiducial galaxy. The line 
profiles are plotted against observed wavelength and also 

Va- 

For both the dynamic and simple models, the IGM out- 
side of the virial radius, though highly ionized out to the 
Stromgren radius, has sufficient optical depth to completely 
suppress the Lyman a line with resonant scattering from the 
Doppler core of the absorption cross-section. Thus the wave- 
length range over which the line is completely suppressed 
corresponds to Va > i;(-Rvir), where u(-Rvir) is the veloc- 
ity of the IGM at the virial radius. For the dynamic model 
v{RviT) = —160 kms~^, in contrast to 15 kms~^ for the 
simple model. 

Because IGM gas is assumed to ionize completely as 
it crosses the virial radius, the line redward of ?;(i?vir) is 
not scattered by the Doppler core of the absorption profile 
of any IGM gas. However, the observed lines are still de- 
pressed compared to the intrinsic profile at Va < v(i?vir) due 
to scattering by the fully neutral IGM outside the Strom- 
gren sphere, via the red damping wing of the absorption 
cross-section. In summary, the optical depth at some wave- 
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Figure 7. Density and ionization state of the IGM as a function 
of radius and velocity, (a) Dynamic IGM model, (b) Simple model. 
The dashed curves are njj, the total (Hl+Hll) hydrogen number 
density; the dotted curves are xjji, the neutral hydrogen fraction; 
and the solid curves are nni = I'Hjnjj. Both plots are for the 
fiducial galaxy with radius in physical Mpc and velocity in the 
frame of the galaxy. 
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Figure 8. Lyman a line profiles. The dotted curve is the intrinsic 
Lyman a line profile of the fiducial galaxy. The solid curve is the 
observed Lyman a line profile using the dynamic IGM model. The 
dashed curve is the observed line using the simple IGM model. 
The velocity where the Lyman a resonance corresponds to the 
wavelength on the bottom axis is shown as Va, on the top axis. 



length (and thus Va) can be thought of as a sum of the 
optical depth due to Doppler core scattering by IGM gas 
at V = Va, and optical depth due to damping wing scat- 
tering by IGM gas at v ^ v{Rs). In Fig. HI the differ- 
ence between the dynamic and simple model predictions at 
Va = —100 kms~^ is primarily a consequence of the presence 
of neutral gas ai v — —100 kms~^ in the dynamic model. At 
Va = —300 kms~^, the difference is entirely a consequence 
of neutral gas between -|-225 and -1-375 kms~^ (the velocity 
range over which the IGM is neutral in the dynamic model 
but ionized in the simple model). 

The integrated observed line flux for the dynamic model 
is only 1.6 per cent of the intrinsic line, compared to 11 per 
cent for the simple model. Thus incorporating a dynamic 
model for the IGM is an important consideration in con- 
necting theoretical predictions for Lyman a emitters into 
observed counts, and for interpreting the ionization state of 
the IGM through Lyman a emission line observations. 



4.3 Galaxy/halo parameter dependence 

In the previous section we showed that our dynamic IGM 
model has a substantial impact on the observed Lyman a 
line from a galaxy at high redshift. In this section we il- 
lustrate how the observed line depends on properties of the 
galaxy and its halo. The parameters we explore are the halo 
mass, galaxy star-formation rate, age of star formation, es- 
cape fraction of ionizing radiation, and relation between the 
velocity width of the emitted line and the halo circular veloc- 
ity. All figures are for the dynamic IGM model only; varying 
the parameters would have analogous effects under the sim- 
ple IGM model. 

For each parameter, we show plots analogous to Fig. Q 
and Fig.|H| as appropriate. From these plots the virial radius 
and Stromgren radius can be read off, as well as the corre- 
sponding IGM infall velocities. The line profile figures show 
the effect of our model parameters, if any, on the intrinsic 
line profile, and the effect on the observed line profile. 

Halo mass In Fig. |U]we recompute the IGM properties for 
the fiducial model, except that we set the halo mass to 10^" 
or 10^^ Mq. Though more massive halos have larger val- 
ues of J?vir (the radii at which the curves terminate), they 
have more negative values of v^R^ir). More massive halos 
also have slightly higher IGM densities around them. The 
combination of larger virial radii (ionizing fiux falls off as 
r~^) and slightly higher IGM densities leads to larger neu- 
tral fractions within the Stromgren spheres of more mas- 
sive halos. The Stromgren spheres of more massive halos 
are slightly smaller (due to larger densities within them), 
but v(Rs) decreases more dramatically with increasing halo 
mass. 

In Fig. 1101 we show the corresponding intrinsic and 
observed Lyman a line profiles. More massive halos have 
broader intrinsic and observed lines (because ud = fvVc and 
Vc increases for more massive halos). A new prediction of 
the dynamic IGM model is that the offset between the in- 
trinsic line center and the blue edge of the observed line (or, 
similarly, the centroid of the observed line) increases with in- 
creasing halo mass. This conclusion is somewhat dependent 
on the assumption that the intrinsic Lyman a emission line 
is centered on the systemic velocity of the galaxy. In practice 
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Figure 9. As in Fig. 13 but note different radial scale. The curves 
are truncated at _Rvir- (a) Halo mass M = 10^" Mq. (b) M = 
10" Mq (fiducial model), (c) M = lO^^ Mq. 
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Figure 11. As in Fig. 13 (a) Star-formation rate tp = 1 Mq yr ^ . 
(b) V = 10 MQyr-i (fiducial model), (c) ^ = 1000 MQjr-^ 
See Fig. I14t for -0 = 10 MQyr"'^. (d) Age of star formation 
tgp = 10'' yr. (e) tsF = 10** yr (fiducial model), (f) tsF = 10^ yr. 
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Figure 10. As in Fig.|S] but note different axis scales. The labels 
show M (in Mq) for the corresponding intrinsic and observed line 
profiles. The fiducial model M = 10^^ Mq is the lighter curve, 
and the intrinsic fiducial line profile is not shown. 



obtaining the intrinsic Lyman a line shape may be impossi- 
ble, but a more complicated model of the intrinsic line may 
be passed through our IGM model to produce a predicted 
observed line profile, which will always have a lower limit in 
wavelength corresponding to Va = «(J?vir). 

Star-formation rate In Figs. 1111 and 1121 we recompute 
the IGM properties and Lyman a line profiles for the 
fiducial model, except that we set the SFR to 1, 100, or 
1000 Mq yr"\ Increasing the SFR increases both the size of 
the Stromgren sphere and the ionization fraction inside. If 
we decrease the SFR from the fiducial value of 10 Mq yr^^ 
to 1 Mq yr~^ then the neutral IGM encroaches all the way 
to v{Rs) = 25 kms~^ (compared to v{Rs) — 225 kms~^ in 
the fiducial case), and thus the damping wing almost com- 



pletely wipes out the observed line. Increasing the SFR to 
100 Mq yr~^, the Stromgren sphere expands out to v{Rs) = 
650 kms~^. Though there is a decrease in the neutral frac- 
tion within the Stromgren sphere, this effect is not suffi- 
cient to lift the complete suppression of the observed line at 
Va > w(7?vir). However, the greater distance (in velocity) to 
the fully neutral IGM means that the observed line redward 
of of Va — v{Rvir) Is less damped than in the fiducial case. 
A further increase of the SFR to 1000 Mq yr~^ expands the 
Stromgren sphere yet further, resulting in very little damped 
absorption at Va < v{Rvir)- Moreover, the neutral fraction 
inside part of the Stromgren sphere is now small enough 
that scattering by the cross-section Doppler core no longer 
completely suppresses the line at Va > w(J?vir). Note the 
notch in the spectrum at Vg = v(Rvii): this is the fe ature 
predicted in QSO spectra bv lBarkana &i Loebl ll2003all . 

Age of star formation In Figs. Illl and ll3l we recompute 
the IGM and line profile properties for the fiducial model, 
except that we set the age of star formation to lO'' or 10^ yr. 
Changing tsF changes the size of the Stromgren sphere, but 
does not change the ionization balance within it. Thus the ef- 
fect on the observed line profile is only to change the amount 
of damping of the line at Va < ii(i?vir). 

Escape fraction of ionizing photons In Fig. 1141 we re- 
compute the IGM properties for the fiducial model, except 
that we set /esc = or 1. Scaling /esc has the same effect 
as scaling the star-formation rate by the same factor, so the 
fiducial model modified to /esc = 1 is the same as the fidu- 
cial model modified to ip = 100 Mq yr~^ (see above). In the 
case of vanishing escape fraction, for any reasonable value of 
the SFR the line will unobservable, because a neutral IGM 
down to v{Rvir) allows only 4 x 10~^ of the intrinsic line to 
be observed. See also Section 
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Figure 12. As in Fig. |S] but different vertical scale. The labels 
show star-formation rate (in Mq yr~^) for the corresponding ob- 
served line profiles. The fiducial model i/> = 10 MQyr~^ is the 
lighter curve, and the intrinsic fiducial line profile is the dotted 
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Figure 13. As in Fig. 13 but different vertical scale. The labels 
show the age of star-formation (in yr) for the corresponding ob- 
served line profiles. The fiducial model tgp = 10* yr is the lighter 
curve, and the intrinsic fiducial line profile is the dotted curve. 
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Figure 14. As in Fig. 13 (a) Ionizing photon escape fraction 
/esc = 0. (b) /esc = 0.1 (fiducial model), (c) /osc = 1. (d) Ionizing 
background J21 = IQ-^. (0) J21 = 1. (f) J21 = 100. 



Velocity width factor In Fig. ll5l we recompute line profile 
properties for the fiducial model, except that we set = 1 
or 2 (/„ does not infiuence the IGM properties). Larger val- 
ues of fv broaden the line, resulting in a greater observable 
fraction of the intrinsic line. In the limit of much larger val- 
ues of fv, more of the line would be transmitted, but the 
line width would be so large that spectroscopic detection of 
the line would become more difficult than for a narrower 
line with less fiux. However, for the values of /„ we consider, 
which should bound the actual values, the observed line is 
still relatively narrow and larger values of /„ produce more 
easily observed lines. 



4.4 IGM ionization dependence 

In this section we vary the ionization state of the IGM to 
show how this infiuences the observed Lyman a line. First we 
will assume, in contrast to the IGM models considered so far, 
that there is a mean ionizing background in the IGM that 
contributes to the ionization of the IGM around the halo, 
in addition to the contribution from the galaxy itself. We 
assume that the universe was reionized before star formation 
turned on in the halo (that is, reionization occurred longer 
than tsF ago), and since then the background has been fixed 
at 



J21 



10- 



-2 Hz 



10« 



(21) 



that 



For the fiducial values oi z — 6.5 and tsF 
corresponds to reionization before z = 7.2. 

In Figs. I14l and ll()l we show the properties of the IGM 
and Lyman a line profiles for values of the ionizing back- 
ground that bracket the physically interesting regimes. First 
we discuss J21 = 10~^, a value so small that the IGM ioniza- 
tion state inside the galaxy's light sphere is completely dom- 
inated by the galaxy's ionizing radiation. The light sphere 
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Figure 15. As in Fig. |S] but different vertical scale. The labels 
show the conversion factor /„ between Doppler line width and 
halo circular velocity for the corresponding intrinsic and observed 
line profiles. The fiducial model = 1.5 is the lighter curve, and 
the intrinsic fiducial line profile is not shown. 
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Figure 16. As in Fig. |S] The labels show ionizing background 
strength J21 for the corresponding observed line profiles. The 
fiducial model with no ionizing background is the lighter curve, 
and the intrinsic fiducial line profile is the dotted curve. 



has a radius of ctsF = 31 pMpc for tsF ~ 10* yr, corre- 
sponding to an IGM velocity of almost 2.5 x 10* kms~^. 
Consequently, the IGM ionization state outside of the light 
sphere is irrelevant. 

Even though the ionizing background does not con- 
tribute at all to the ionization rate (see eq. 118^ at any ra- 
dius of interest, the pre-reionization of the universe before 
the galaxy turned on allows its Stromgren sphere to expand 
at near light speed. Consequently even the IGM at the ve- 
locity of 1200 kms^^ has a neutral fraction of 0.01. By the 
Gunn-Peterson trough test, this IGM wouldn't even qual- 
ify as reionized. However, this level of ionization completely 
eliminates any damping wing optical depth. Thus the ob- 
served line profile follows the intrinsic line profile from the 
red wing up to Va = v{Rvir), then goes to zero as in the 
fiducial case. This is the observed behavior for any ioniz- 
ing background w ith J21 < 1, the range observed at 2; < 5 
JScott et al.l l2000^. 

Ionizing backgrounds in the range J21 ~ 10-100 
strongly ionize the mean density IGM. Since the IGM den- 
sity decreases from a local maximum at iivir toward the cos- 
mic mean as a function of increasing radius, IGM gas just 
outside of the virial radius still strongly scatters the intrinsic 
line at Va corresponding to the velocity of the overdense gas. 
As the density decreases, though, the neutral fraction drops 
as well, and the combination leads to greater transmission of 
the intrinsic line toward the blue side of the line. This would 
generate an unusually-shaped observed line with two peaks, 
one of which lies to the blue of the intrinsic line center. 

At yet higher values of J21 > 10"^, the IGM is so strongly 
ionized that the entire intrinsic line would be observed. We 
note again that we expect J21 < 1, though we are unaware 



of any predictions for J21 at z > 6 in models that reproduce 
the WMAP optical depth value. 

Alternately, rather than considering a fixed ionizing 
background in a reionized universe, we can fix the neu- 
tral fraction of the IGM outside of the galaxy's Stromgren 
sphere, a^tn'^. This approach is more useful for considering 
a weakly-ionized IGM with x^^^^ ~ 0.5 ^e.g.. ICenI EoO.^al . 
where a fully reionized universe almost entirely recombines) . 
Figure [T7I shows the observed Lyman a line profiles of the 
fiducial galaxy when is set to 0.5 and 0.05. As ex- 

pected, the decreased Hi abundance in the 'neutral' IGM 
means that the red wing of the Lyman a line is less sup- 
pressed by damping wing scattering. There is no effect on 
the Doppler scattering that scatters the center of the intrin- 
sic line. See also Section f6. 21 

4.5 Redshift dependence 

Last, we turn to the effects due solely to changes in redshift. 
In the future there will (hopefully) be surveys for strong Ly- 
man a emission lines at several different redshifts at z > 6. 
We study redshift dependence by keeping the other free pa- 
rameters fixed; realistically, other parameters, such as star- 
formation rate, age of star formation, escape fraction, and 
the typical halo mass, may also depend on redshift. Compar- 
ison of the galaxy/halo parameter dependencies described 
above in conjunction with the effects shown here can be used 
to predict trends for a specific model of the evolution of the 
galaxy/halo parameters. Disentangling the redshift depen- 
dence of the IGM ionization state from the other redshift 
dependences may still be very difficult in practice. 

Proceeding with simple cases to understand the ex- 
pected effects, we place the fiducial halo at 2 = 8.8 and 
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Figure 17. As in Fig. |S] but with different vertical scale. The 
labels show the neutral fraction of the IGM outside of the galaxy's 
Stromgren sphere for the corresponding observed line profiles. 
The fiducial model with a neutral IGM is the lighter curve, and 
the intrinsic fiducial line profile is the dotted curve. 
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Figure 18. As in Fig. |7| but note different radial and velocity 
scales. The curves are truncated at i?vir- (a) Redshift z = 6.5 
(fiducial model), (b) z = 8.8. (c) z = 17. 



z = 17. Figures 1181 and 1191 show the IGM and Lyman a 
line profile properties. Because we keep the mass of the halo 
fixed, the virial radius shrinks at increasing redshift. This 
results in an increasingly negative velocity of the IGM at 
the virial radius. The IGM is also denser at higher redshift, 
resulting in higher neutral fractions when the ionization bal- 
ance is calculated. An additional effect is the increased cir- 
cular velocity of the halo, which broadens the intrinsic line. 
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Figure 19. Lyman a line profiles. The wavelength scale is ex- 
pressed as an offset from the Lyman a resonance in the frame 
of the galaxy, chosen so that the scale is the same as in Fig. 
ISl The labels show the redshift for the corresponding observed 
line profiles. The narrower intrinsic profile (dotted curve) is for 
z = 8.8; the wider intrinsic profile (dotted curve) is for z = 17. 
The fiducial model z = 6.5 is the lighter curve, and the intrinsic 
fiducial line profile is not shown. 



far. High velocity (~ 360 kms~^) redshifts of Lyman a emis- 
sion compared to other inter stellar features are o bserved in 
z = 3 Lyman Break Galaxies iShaplev et al .120031) : QSO ab- 
sorption line cross-correlation suggests that these winds may 
have blown to large distances (~ 0.2 pMpc, Adclbcrgcr ct ajj 



1200: 

Galactic winds, if present in galaxies at 2; > 6, may 
have two important consequences for the interpretations of 
Lyman a emission from galaxies during reionization. The 
first is that as the wind blows through the IGM, it may col- 
lisionally ionize the hydro gen it passes, increasin g the ion- 
ization fraction of this gas llAdelberger et al.l2003l) . The sec- 
ond effect is that the intrinsic Lyman a line may not have 
a Doppler shape centered on the galaxy's redshift. Instead, 
it may already display an asymmetric profile centered con- 
sider ably to the red of the systemic redshift iShaplev et alJ 
I2OO3I where they find a mean redshift of 360 kms~^). In- 
spection of Fig. |S| shows that the IGM optical depth for 
the fiducial galaxy is only ~ 1 at redshifts of more than 
200 kms"^"*. 

Since we lack a good physical model to calculate wind 
properties, we adopt a very simple assumption to illustrate 
the possible effects of large-scale winds on the IGM (the first 
effect described above) . We assume that a wind blows out to 



5 GALACTIC WINDS 

We now consider another, possibly crucial, alteration to the 
model of Lyman a emitter plus IGM scattering described so 



* Note that in observations velocity (and redshift) is measured 
in the observer's frame, whereas in our plots Va is in the frame of 
the emitting galaxy, hence a line redshifted by 360 kms~^ falls 
at Va = —360 kms"'^. 
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Figure 20. As in Fig. |S] but different vertical scale. The labels 
show star-formation timescale (in yr) for the corresponding ob- 
served line profiles. For each igp, the light solid curve is the profile 
without assuming any winds (as in Fig. 1131 . and the long-dashed 
curve is the observed line profile assuming a wind blows to the 
Stromgren radius. The intrinsic fiducial line profile is the dotted 
curve. 



the Stromgren radius of the galaxy (see Fig. lllH . completely 
ionizing all of the gas within that radius. Figure l2?)l shows 
the observed Lyman a line profiles including such winds. 
The result for the fiducial model is that the observed line 
becomes roughly symmetric, though still centered to the red 
of the intrinsic line center, and the observed flux doubles. 
If tsF is changed to 10^ yr, the Stromgren sphere is much 
larger (see Fig. lllH and the observed fiux increases by almost 
a factor of 5. 

The required wind velocity can be calculated by divid- 
ing the Stromgren radii (illustrated in previous plots) by tsF- 
For the fiducial halo, the Stromgren radius is 0.38 pMpc, 
giving a wind velocity over 10** yr of 3700 kms"\ Almost 
any wind model would predict a smaller wind velocity than 
this, so though Fig. I20l shows the qualitative effect, even the 
10* yr curve is an overestimate of the true effect, assum- 



ing g alactic winds blow to 0.2 pMpc, as in lAdelberger et alJ 
ll2003 ). 

Winds don't affect the damping wing suppression of 
the line from absorption by neutral (in the case of a non- 
reionized IGM) gas outside of the Stromgren sphere (the 
solid and long-dashed curves match at Vc, < v{Rvii)). Thus 
the depression of the observed line profiles compared to the 
intrinsic line profiles of plots like Fig.|S|can be used to esti- 
mate the transmission of a highly-redshifted Lyman a line. 
Since the damping wing optical depth is set by the IGM 
outside of the Stromgren sphere, parameters that affect the 
size of the Stromgren sphere (e.g., tsF, Fig. llH . and the ion- 
ization state of the IGM (e.g., J21, Fig. 1141 and i]^*^), have 
a large influence on the transmission at Vc, < w(/?vir). In 
the next section we illustrate the effect galactic winds may 



have on the observed Lyman a line for all of the parameter 
variations discussed in Section 0| 



6 DISCUSSION 

6.1 Lyman a line fluxes 

First we present three plots to summarize the results of the 
previous sections. Figure 1211 shows the intrinsic Lyman a 
line flux (calculated as described in Section 13.11 as open 
triangles for each model listed on the bottom axis. The solid 
circles show the integrated flux of the observed Lyman a line 
calculated from the intrinsic Lyman a line flux and observed 
line profiles presented in Section 2] 

Figure I2H shows the intrinsic Lyman a line fiux as open 
triangles for each model, as in Fig.HH Here the open squares 
are the integrated fiux of the observed Lyman a line calcu- 
lated assuming a galactic wind completely ionizes the IGM 
out to the Stromgren radius of the galaxy (see Section |S| 
this is probably an overestimate of the expected effect, but 
comparison with Fig. 1211 shows that for many models this 
is a small effect regardless). The solid diamonds are the in- 
tegrated flux of the observed Lyman a line calculated as- 
suming the intrinsic Lyman a line is emitted entirely at 
360 kms~^ redward of the systemic Lyman a wavelength 
(see Section |^. The consequence of moving the intrinsic 
line to the red is that the observed line flux is generally 
at least 20 per cent of the intrinsic line flux for most of the 
models (independent of whether of not the wind blows into 
the IGM) . This is illustrated more explicitly in Figure 1231 
However, more massive galaxies with ii(i?vir) ^ 360 kms~^ 
(such as our M = 10^^ M© model) would require Lyman a 
lines even more redshifted than 360 kms~^ to make a dif- 
ference compared to the no-wind model. There is no solid 
diamond for M = 10^^ Mq in either Fig. ^ or Fig. ^ 
because the IGM transmission at Va = —360 kms~^ is for- 
mally zero. A realistic line profile (not concentrated entirely 
at Va = —360 kms~^) would result in a flux similar to the 
no-wind model. 

Figure 1^ shows the transmission of the Lyman a line, 
Ta , deflned as the ratio of the observed to intrinsic Lyman a 
line flux. The shapes of the symbols match their shapes in 
Figs. I2l| and 1221 Quite generally, an intrinsic Lyman a line 
redshifted by 360 kms~^ produces an observed line sup- 
pressed by less than a factor of 5, for galaxies at 2: = 6.5 in a 
neutral IGM. If the IGM is reasonably ionized, a redshifted 
Lyman a line may survive almost unscattered (e.g., solid di- 
amond for ajjji^ = 0.05 model), whereas if the intrinsic line 
is at the systemic redshift in those models, the observed line 
is almost an order of magnitude weaker (e.g., open circle for 



x^S^ = 0.05 model). 



6.2 Observed 2: = 6.5 Lyman a emitters 

To compare observed Lyman a line fluxes to our predictions, 
we need an estimate of the intrinsic line flux. The observed 
Lyman a hne flux provides a lower limit to the intrinsic flux, 
but current observations of 2: = 6.5 do provide one avenue to 
estimate the intrinsic flux. The UV continuum t o the red of 
the Lyman a is a star-formation rate indicator (iKennicutii 
1I9981I : near-IR photometry of z = 6.5 galaxies provides an 
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Figure 21. Intrinsic and observed Lyman a line flux as a function 
of model parameters. Each model plotted is described by a label 
on the bottom axis that refers to the fiducial model in the case 
of 'fiduc,' and refers to the parameter varied from the fiducial 
model in the rest of the cases (see Section Numbers in square 
braces are the logj^g of the parameter value. The open triangles 
are the intrinsic Lyman a line fliax, and the solid circles are the 
observed Lyman a line flux, based on the line proflles calculated 
in Section |1] 



estimate of their rest-frame UV continua. Two of the three 



confirmed z = 6.5 galaxies have SFRs esti mated in this wa: 
10 Moyr"^ for HCM 6A (hereafter 'H.' IHu et all 12002 



i, 



and 36 yr~^ for SDF J132415. 7+273058 (hereafter 'K,' 
iKodaira et al.l 1200311 (in our cosmology). As discussed in 
Section 13.11 converting the SFR to an intrinsic Lyman a 
line flux depends on many parameters, but for our fiducial 
choices we estimate intrinsic Lyman a line luminosities of 
3.4 X 10^^ erg s~^ fo r H and 1.2 x 10'*^ erg s~^ for K Note 
that both lHu et"ai] (|20b2a) and iKodaira et alJ i2003l) esti- 
mate Lyman a lower intrinsic luminosities by us ing an em- 
pirica l relation between Ha luminosity and SFR iKennicutti 
Il998^ and converting Hq luminosity to Lyman a luminos- 
ity. The corresponding intrinsic Lyman a line fluxes are 
6.9x 10"" erg s"^ cm'^ for H and 2.4x 10"^*^ erg s"^ cm'^ 
for K. Comparison with the observed Lyman a fluxes (cor- 
rected for lensing in the case of 'H') gives estimated observed 
to intrinsic ratios of Tq = 0.09 for both H and K. 

From Fig. 1231 we see that appropriate models with a 
neutral IGM and some velocity offset between the intrinsic 
Lyman a line and the systemic galaxy redshift can produce 
values of Tq ~ 0.1. Alternately models with a partially or 
mostly ionized IGM without any wind effect also have Ta ~ 
0.1. 

The redshift of the intrinsic Lyman a line with respect 
to the systemic redshift of the galaxy is a clearly key param- 
eter for interpreting observed Lyman a lines. Fortunately it 
may be estimated from non-resonant recombination lines, 
such as Ha or H/3. For galaxies at z > 6.5 these lines fall 
at A > 3.6 ^m, unobservable from the ground in the near- 
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Figure 22. Intrinsic and observed Lyman a line flux as a func- 
tion of model parameters, including possible effects due to galactic 
winds. The model labelling is described in Fig. 1211 The open tri- 
angles are the intrinsic Lyman a line flux, as in Fig. 1211 The open 
squares are the integrated flux of the observed Lyman a line cal- 
culated assuming a galactic wind completely ionizes the IGM out 
to the Stromgren radius of the galaxy. The solid diamonds are the 
integrated flux of the observed Lyman a line calculated assum- 
ing the intrinsic Lyman a line is emitted entirely at 360 kms~^ 
redward of the systemic Lyman a wavelength. See Section l5l 



future. The James Webb Space Telescope [JWST) is pro- 
jected to have sufficient sensitivity to detect Ha and H/3 
from the z = 6.5 galaxies discovered so far. Measurements 
of Ha and H/3 will also provide another SFR estimate (with 
/osc no longer entering into the conversion between SFR and 
intrinsic Lyman a line flux), as well as a dust extinction es- 
timate. 

An example of an optimistic future observation would 
be a measurement of Ha and H/3 line redshifts and fluxes 
for a z = 6.5 galaxy with a measured Lyman a line and 
UV continuum. If the ratio of Ha to H/3 line flux matches 
the prediction for no extinction, and the SFR deduced from 
from those Balmer lines agrees with the SFR estimated from 
the UV continuum, then a reasonable conclusion would be 
that dust extinction is not very important in the galaxy^. 
Proceeding under the assumption of a normal IMF and no 
extinction, the intrinsic Lyman a line flux may then be esti- 
mated. If the Lyman a line were found to be centered near 
the Balmer line redshift (assumed to be the systemic redshift 
of the galaxy), then we could compare the measured value 
of Ta with the circles in Fig. 1231 at that point a Ta could 
be converted into an estimate of the neutral fraction of the 
IGM, though with plenty of uncertainty associated with un- 
known parameters such as M, tsF, and /„. A more clear-cut 



^ This conclusion could be strengthened by a stringent upper 
limit on, e.g.. Hen 1640 A, a line that would be present for an 
IMF weighted toward very massive stars, (^Schaorcr 2003). 
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Figure 23. Ratio of observed to intrinsic Lyman a line flux as a 
function of model parameters, for botli no wind and wind models. 
The model labelling is described in Fig. 1211 The point shapes 
are the same as in Figs. En and El open circles show no-wind 
models, open squares show models where a strong wind clears the 
Stromgren sphere of IGM, and solid diamonds show models where 
the intrinsic Lyman a is concentrated at 360 kms^^ red ward of 
the systemic redshift. 



conclusion may be reached if the Lyman a line is found to 
be substantially redshifted compared to the Balmer redshift 
and the measured Ta is close to 1. This would be a strong 
signal that the IGM was mostly ionized at that redshift (see 
below), and may moreover place an interesting limit on the 
mass of the halo, depending on the shape of the observed 
Lyman a line. 

In advance of JWST, Fig. l23l makes it clear that drawing 
any firm conclusions regarding the ionization fraction of the 
IGM based only on the observed Lyman a lines and UV 
continua of z = 6.5 is very difficult. 

Before moving on, we present two more plots of interest 
regarding the observability of Lyman a lines. We have seen 
that large Stromgren spheres around galaxies allow more of 
the intrinsic Lyman a line to be observed (e.g.. Fig. 
The ionization of the IGM is directly proportional to the 
parameter /esc, so the largest possible Stromgren sphere is 
created when /esc = 1- However, since the intrinsic Lyman a 
line flux is proportional to (1 — /esc), the strongest intrin- 
sic Lyman a line is produced when /esc = 0. Clearly the 
strongest observed line comes from some intermediate value 
of /esc- In Fig. 1241 we plot the intrinsic and observed Ly- 
man a line flux for two models, the fiducial model, and the 
t/j — 100 Mq yr~^ model. The observed line flux in both 
models peaks near /esc — 0.35. 

Once reliable estimates of Ta are available, they can be 
used to estimate the neutral fraction of the IGM in the in- 
teresting regime of xui ~ 0.5. Figure 1^ illustrates as a 
function of the model parameter (see Section [4.411 . for 

the fiducial galaxy. The solid curve was calculated assum- 



01 




Figure 24. Intrinsic and observed Lyman a line flux as a func- 
tion of ionizing photon escape fraction. The dotted curves are 
the intrinsic Lyman a line flux; the solid curves are the observed 
line fluxes. The lower curves are for the flducial model with /esc 
varied; the upper curves are for the flducial model modified to 
= 100 M0yr-l. 



ing a Doppler profile intrinsic Lyman a line centered at the 
galaxy's systemic redshift, and the long-dashed curve was 
calculated assuming the intrinsic Lyman a line is concen- 
trated 360 kms~^ redward of the systemic velocity. While 
an accurate measurement of Ta can constrain the IGM neu- 
tral fraction, this is only possible with reasonable knowledge 
of the redshift of the intrinsic Lyman a line compared to the 
systemic velocity of the galaxy. 



6.3 Model assumptions 

Now we return to a discussion of the assumptions made in 
this paper and their possible consequences on our conclu- 
sions. First we discuss the possible effects of dust in a Ly- 
man a emitting galaxy. If dust absorbs ionizing photons, 
then the star-formation rate in our models can simply be 
rescaled to take that into account. Alternately, dust may ex- 
tinct the Lyman a line and UV continuum almost equally, 
but extinct the Balmer lines considerably less. In this case 
the Lyman a line extinction can be estimated and corrected 
for. Dust may also affect the UV continuum differently from 
the recombination emission, in the event that dust is dis- 
tributed differently with respect to stars and H ll regions. In 
this case the extinction value measured by the Balmer lines 
alone (with an extinction curve) may be used to correct the 
Lyman a line, assuming that resonant photons are extincted 
no differently than non-resonant photons. 

Dust extinction can cause great ambiguity in the inter- 
pretation of Lyman a lines, even when UV continuum and 
Balmer lines are measured, because resonant Lyman a pho- 
tons may be extincted differently from non-resonant UV con- 
tinuum and Balmer line photons. Dust may extinct Lyman a 
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Figure 25. Ratio of observed to intrinsic Lyman a line flux as 
a function of a;]^'^. The solid curve assumes the intrinsic line 
profile described in Section 13.11 The long-dashed curve assumes 
the intrinsic Lyman a line is concentrated 360 kms~^ redward 
of the systemic velocity (see Section |FJ. Both curves are for the 
fiducial galaxy. 



photons preferentially over UV continuum photons, due to 
the extra path-length through the galaxy that resonantly 
scattered photons travel before escaping. Certain gas and 
dust g eometries, how ever, may mitigate or even reverse that 
effect jNeufelJl99lh . There is no direct observable to mea- 
sure this possible effect, except that an extinction-corrected 
SFR estimated from Balmer lines does place an upper limit 
on the intrinsic Lyman a line flux. 

More fundamental assumptions are associated with our 
model of the IGM around Lyman a emitting galaxies. We as- 
sumed the presence of an accretion shock a t the virial radius 
of the halo containing the galaxy (see also lBarkana fc Loebl 
l2003d) . If these galaxies hav e no a ccretion shock, as sug- 
gested by iBirnboim &: Dekell i2003ft , then our conclusions 
are greatly changed. It is the presence of the accretion shock 
that generates the characteristic wavelength scale where line 
transmission begins (see Section 14.211 . In the absence of a 
virial shock, the infalling gas will eventually shock when it 
reaches the galaxy. The observed Lyman a line would be 
suppressed blueward of that velocity, implying only intrin- 
sic Lyman a lines that were very redshifted with respect 
to the systemic redshift of the galaxy could be observed. 
Another consequence of the lack of a virial shock is that Ly- 
man a emission from the co oling of the IGM gas wo uld be 
located close to the galaxy iBirnboim fc DekeJl2003|) . thus 
any observed Lyman a line may have a contribution from 
contribution from cooling radiation as well as star-formation 
powered Lyman a. If future observations show an extremely 
redshifted Lyman a lines at high redshifts, this may be in- 
dicative of cooling radiation from IGM gas falling unshocked 
down to the center of a halo. 

In Section EHU we discussed our assumptions about 



the density distribution of the IGM around the Lyman a 
emitting galaxy. We calculate the ensemble-averaged IGM 
density distribution; even for the spherically-averaged den- 
sity profile, deviations from this average are expected from 
galaxy to galaxy. Additionally, we estimated that structure 
near the emitting galaxy has only a small effect on the den- 
sity distribution of the IGM. However, we did not address 
the effect of collapsed and collapsing structures on the ve- 
locity field of the IGM, which is crucial (more so than IGM 
density) to our predicted line proflles. Numerical simula- 
tions may be used to estimate both the variance in the IGM 
structure from galaxy to galaxy, and also the density and ve- 
locity distributions of the IGM along the line of sight toward 
a galaxy accounting for substructure. 



7 SUMMARY 

We have calculated observed Lyman a line profiles for high- 
redshift galaxies, with a focus on the prospects of using ob- 
served Lyman a lines to probe the ionization of the universe 
at z > 6.5. Toward this end we have investigated many of 
the parameters that influence the observed Lyman ct line 
from a high-redshift galaxy. We conclude that future mea- 
surements with JWST may permit the estimation of IGM 
neutral fractions aX z > 6.5, but current z — 6.5 data do 
not yet place a strong constraint on the neutral fraction at 
that redshift if we allow consideration of galactic winds. If 
we take at face value our estimates of the integrated Ly- 
man a transmission of the two z = 6.5 galaxies for which 
estimation is possible, and assume our simplest model (with 
no winds), we would conclude that the IGM neutral fraction 
is xm ^ 0.1 at z = 6.5. 

Interpretation of an observed Lyman a line alone (or 
even with a measurement of the UV continuum) is extremely 
difficult, due to the host of parameters that influence it. 
However there is hope that with the rest-frame optical spec- 
tra that JWST will provide for high-redshift galaxies, Ly- 
man a emission lines will yet prove fruitful for constraining 
the reionization history of the universe at z > 6.5 
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